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Abstract 

Line-by-line calculations are becoming the standard 
procedure for carrying spectral simulations. However, it is 
important to insure the accuracy of such spectral simula- 
tions through the choice of adapted models for the simula- 
tion of key parameters such as line position, intensity, and 
shape. Moreover, it is necessary to rely on accurate spec- 
tral data to guaranty the accuracy of the simulated spec- 
tra. A discussion on the most accurate models available 
for such calculations is presented for diatomic and linear 
polyatomic discrete radiation, and possible reductions on 
the number of calculated lines are discussed in order to re- 
duce memory and computational overheads. Examples of 
different approaches for the simulation of experimentally 
determined low-pressure molecular spectra are presented. 
The accuracy of different simulation approaches is dis- 
cussed and it is verified that a careful choice of applied 
computational models and spectroscopic datasets yields 
precise approximations of the measured spectra. 
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1 Introduction 

Advances in the capabilities of modern day computing 
have allowed carrying line-by-line calculations in a sys- 
tematic manner for the field of spectroscopy modelling. 
Various applications such as species concentrations and 
temperature measurements in low pressure plasmas 012], 
or the calculation of radiative fluxes for atmospheric en- 
try flows 0, greatly benefit from such techniques. How- 
ever, the accuracy of such spectral simulations may vary 



consequently, depending on the methods used for line-by- 
line calculations, but also on the applied spectroscopic 
datasets. 

The different methods used for the calculation of key 
spectral parameters such as line positions, intensities, and 
shapes are discussed in this paper. Available state of the 
art models for each of these parameters are presented, 
and some simplifications to such spectral models are dis- 
cussed, leading to lower memory requirements and calcu- 
lation times for the used computing systems. The mod- 
els presented and discussed in this paper are valid for 
diatomic rovibronic transitions, but also for linear poly- 
atomic rovibrational transitions (such as transitions from 
the CO2 molecule). The second part of this work presents 
some simulations of experimentally measured high res- 
olution spectra issued from low pressure and high en- 
thalpy plasmas. The discrepancies that may derive from 
the selection of different simulation models and spectral 
datasets will be highlighted through the comparison of 
different simulated spectra with the measured spectra. It 
will be verified that a careful selection of adequate mod- 
els, linked to a selection of accurate spectroscopic data, 
may yield a very accurate reproduction of high resolution 
measured spectra. 

2 Theoretical Models for the Simulation 
of Bound Spectra in the Line-by-Line 
Approach 

Discrete molecular radiation can be characterized un- 
ambiguously through three parameters: line position, in- 
tensity, and shape 

Determination of line positions depends on the quan- 
tification of the energy levels of the molecule bound states. 
Line intensities depend on the probabilities of transition 
between the different states, as well as on the population 
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of these states. Finally, line shapes depend on the local 
conditions of the gas in which the transition takes place 1 . 

An overview of the different methods and approaches to 
several levels of accuracy will be presented in this section 
for these three parameters. 



exist, and care should be exercised when using published 
coefficients. 

Traditional spectroscopic developments are related to 
Zare's Klein-Dunham coefficients by: 



2.1 Line Positions Calculations 

Three different approaches exist for the calculation of 
quantum level energies which will determine line positions 
of bound molecular spectra: 

The first one allows a broad calculation of any number 
of vibrational and rotational levels for a given electronic 
transition of a molecule. Calculations are performed us- 
ing equilibrium constants, which give the vibrational and 
rotational constants for the given electronic transition in 
a broad range of vibration and rotation levels. 

The second one uses level-by-level spectroscopic con- 
stants and allows the calculation of any number of rota- 
tional levels for given electronic and vibrational transition 
levels of a molecule, as the band origin and rotational con- 
stants are set for each vibration level. This approach gen- 
erally allows a better determination of specific vibration 
levels energies but prevents one from simulating further 
levels than those for which spectroscopic data is available. 

The third one, and the more accurate, requires diag- 
onalising the corresponding hamiltonian matrix for each 
rovibronic state. This approach can be useful when the 
experimental spectra to be simulated is strongly per- 
turbed although it leads to larger computational times. 

2.1.1 Line position calculations using equilibrium 
constants in matrix form 

Klein-Dunham coefficients allow a clear and unambigu- 
ous determination of level positions, compared to the tra- 
ditional spectroscopic developments, prone to confusions 
and errors (for instance, the parameter j e , useful for the 
calculation of the rotational constant B v can sometimes 
be confused with the spin-rotation interaction coefficient 
7). 

Level positions (in cm 1 ) are calculated according to 
the following relation: 



E e , v ,j = Yu (v + 1/2)* [F( J)] j 



(1) 



The formalism proposed here is consistent with the ef- 
fective Hamiltonian of Zare Other formalisms may 

1 line positions shifts can also occur from broadening mech- 



G v = F i0 (« + i) =T e +cu e (v + ^) (2a) 



B v = J2 Yl ( v+ ^) =B e -a e (v + ^) (2b) 



+ r Ye\V + 



D v = - Y ^ ( v +1) =D e + pJv+^)+. 



(2c) 
(2d) 



The use of Klein-Dunham expansions, unlike tradi- 
tional developments, may also prevent situations were 
neglecting higher order corrections leads to considerable 
shifts for the position of calculated lines compared to the 
experimental spectrum 0. 

Setting a matrix of Klein-Dunham coefficients such as 
i = 10 and j = 7 in computer routines for line position 
calculations suffices in order to account for the available 
polynomial expansions as verified by a broad review of 
available spectroscopic coefficients by the author. 

For multiplet transitions, the expression for the level 
energies differs slightly from Eq. Q but the general form 
of the Klein-Dunham matrix can be used. Expressions 
for the different multiplet level energies are presented in 
appendix lAl 

2.1.2 Further line position calculation methods 

When the level positions can no longer be accu- 
rately approximated through Klein-Dunham expansions 
(as when vibrational perturbations of the spectra are 
present), spectroscopic constants for each vibrational 
level must be used. Level spectroscopic constants ob- 
tained from fits of the rotational lines for each vibrational 
band are given in this case. Using such level constants 
usually results in more accurate predictions of the level 
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energies. However calculations are restricted to the vibra- 
tional levels were spectroscopic constants are available, 
unlike Klein-Dunham expansions which allow higher level 
extrapolations of the available spectroscopic data. 

When perturbations are present in the spectra, the 
polynomial expansions described previously no longer suf- 
fice for the accurate simulation of line positions. Instead, 
one has to solve the Hamiltonian matrix, taking into ac- 
count the effects of the perturbing states using the pertur- 
bation method This leads to very precise calculations 
of the line positions (typically less than 0.1 cm -1 ). How- 
ever, this method requires the calculation of the proper 
values of a n * n matrix for each rovibronic state where n 
is the level multiplicity. 

2.2 Line Intensities Calculations 

An overview of the different methods used for the cal- 
culation of line emission and absorption as well as the 
different difficulties and approximations seldom encoun- 
tered will be presented in this section. 

The emission coefficient for a single line (excluding 
broadening mechanisms) is calculated according to the 
relation 



N u A ut AE 

47T 



J' 



(•3) 



Eq. |31 highlights the additional difficulties related to 
line intensity calculations. These depend not only on the 
line positions (through the accounting of the transition 
energy A_E), but also on the transition probabilities A u i 
and the number density of the initial state N u . Although 
the two former quantities only depend on the transition 
parameters, being calculated according to quantum me- 
chanics laws, the latter depends on the state of the studied 



gas. 



2.2.1 



Dependence between line emission and 
absorption coefficients 



The radiative properties of a gas can be unambiguously 
known through the determination of it's wavelength- 
dependent emission and absorption coefficients e„ and 
a(v). These two quantities are not independent however, 
and line absorption coefficients can be determined from 
the line emission coefficients. 



The absorption coefficient including spontaneous and 
induced absorption (adopting the normalization factor 
2J + 1 for the Honl-London coefficients and excluding 
broadening mechanisms) is given by the relation 



hi' 



a(v) = — {NiBiu — N u B u i) 
c 



(4) 



In complete thermodynamic equilibrium (CTE), the 
emission and absorption coefficients are related through 
Planck's law: 



^ =Ll=2hcu3 rwr 1 \ (5) 

One can therefore deduce the following relationships 
between the Einstein spontaneous emission coefficient 
A u i, the Einstein spontaneous absorption coefficient Bi u , 
and the Einstein induced absorption coefficients B u f. 



B u i 



and 



giBiu — guB u i 



As these relationships hold even in non-local thermody- 
namic (NLTE) conditions, after simple algebraic manip- 
ulation, one can relate the line emission and absorption 
coefficients (in wavenumber units) through the expression 



ot(v) = Bp- 



' 2hc 2 v 3 \g u n u ) 
which means that only one of these two coefficients 
needs to be calculated. 

For the more restrictive case of a Boltzmann dis- 
tribution of the internal levels populations where 
Tit = gi exp ( — &) IQtot one recovers the relation be- 
tween the emission and absorption coefficient defined by 
Planck's Law (Eq. 0. 

2.2.2 Determination of the initial quantum levels 
populations N u 

In thermodynamic equilibrium, the population of a 
quantum level can be calculated straightforwardly from 
the well-known Boltzmann equilibrium relation 



Nt 
N 



(7) 



with Qi = gi exp ( - 



However, accurate calculations of atomic and molecular 
species partition functions requires a set of accurate spec- 
troscopic constants up to higher quantum levels (close to 
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the dissociation limits for molecules and ionization limits 
for atoms) . Namely, the lowering of the ionisation thresh- 
old for atomic species has to be considered in the plasma 
state, when including the contribution of the atomic Ry- 
dberg states to the overall partition function Q. Also, 
for molecular partition functions calculations, it is neces- 
sary to determine the maximum rovibronic levels which 
can be achieved before the dissociation of the molecule 
occurs (including superdissociative states). 

Out of thermodynamic equilibrium, no straightforward 
method for calculating the levels populations exists, and 
one has to resort to state-to-state models which explicitly 
take into account the different possible discrete states of 
a gas species. The development of accurate state-to-state 
models has been carried by different research teams, and 
the reader should refer to the references [S] HO ^] ^3 H2I 
1131 1111 for a more detailed description of such models. 

2.2.3 Determination of the transition probabilities A u i 

As it was seen before, the global transition probability 
can be decomposed into a vibronic and rotational part: 



A u l — A%n v ll ■ A%fjll 



(8) 



The rotational transition probability is calculated ac- 
cording to the Honl-London factors which depend on the 
electronic transition type ( n A ^ n A). 



A 



A J & A" J" 



A" J" — 



2,1' + 1 

in which the normalisation rule reads 



(9) 



E S a'< J ;,(J') = (2J' + 1) 



(10) 



Some authors (151 use a slightly normalisation rule 



£ Sa- J ;,(J') = (2 - <5 ,A')(2S + 1)(2J' + 1) (11) 



Expressions of the Honl-London factors for the different 
types transitions are given by several authors and can be 
found in Refs. IT71 ITS1 IT51 15^1 . 

The vibronic part A%i? v ti of the transition probability 
can be written as a function of it's vibronic transition 
moment. This can be written in atomic units as 



= ^rt J 2 ~ do ? +A : ] (k' v "T (12) 



3hc 3 " " 
2.026 • 10~ D z7 



(2- Jo,A') 
6-3 ( 2 — Va'+A") 

(2-<5 ,A') 



As the vibronic transition moment \lft" P v ) cannot be 



resolved for each multiplet transition, a mean value of the 
transition I 3?^ v I moment is rather used: 



(2-5 0iA '+A")(2<? + l) ( 13 ) 

The vibronic transition moment can be calculated for 
each vibronic transition as a function of the electronic 
transition moment and the transition initial and final 
wavefunctions 

($€ v ") 2 = (J VvW(r)VvKr)*V (14) 

Expressions for the electronic transition moment can 
be found in the literature, either from spectroscopic mea- 
surements, or from "ab-initio" calculations. This last 
method is usually preferred, as nowadays, quantum meth- 
ods have achieved a very good precision |2*T) . The vibra- 
tional wavefunctions are determined by solving the radial 
Schrodinger equation on the potential curves of the upper 
and lower level potential curves. 

Potential curves can be either calculated using "ab- 
initio" methods, or reconstructed through the Rydberg- 
Klein-Rees (RKR) HO IH method according to ex- 
perimental spectroscopic data. As modern spectroscopy 
is able to resolve line positions to less than the cm -1 , 
level energies can be known to a greater accuracy than 
using "ab-initio" methods. However, "ab-initio" methods 
are able to reproduce the entire potential curve, whereas 
the RKR method can only yield the region of the poten- 
tial curve where measured data is available. To overcome 
this problem, the central part of the potential curve cal- 
culated through the RKR method is extrapolated by a re- 
pulsive potential at narrower internuclear distances, and 
by a Hulburth and Hirschfelder 1251 potential at larger 
internuclear distances, provided that the state dissocia- 
tion energy is know. This method can only be applied 
for electronic states with a single potential with a shape 
close to a well, but this is fortunately the case for most 
of the electronic states for the molecules encountered in 
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gas spectroscopy. A more detailed overview of the cal- 
culation of potential curves and vibronic wavefunctions 
can be found in 1261 . An example of a calculation of po- 
tential curves (using the RKR method) and vibrational 
wavefunctions (solving the radial Schrodinger equation) 
is presented in Fig. Q 



.x 10 




1.5 2 2.5 

Internuclear Distance (A) 

Fig. 1 : Example of an RKR calculation of the potential curves 
of the upper and lower electronic states for the CN 
Violet transition. The wavefunctions of different vi- 
brational levels of both electronic states are also re- 
produced at an arbitrary scale 



2.3 Simulation of Linear Polyatomic 
Rovibrational Spectra Using the 
Line-by-Line Approach 

Calculations for linear polyatomic species present 
further difficulties compared to diatomic species, as 
molecules have several vibration modes (bending, sym- 
metric/asymmetric stretch, etc.). Emission spectra from 
these molecules results mainly from rovibrational transi- 
tions, and a set of vibrational equilibrium constants can 
still be defined, with a more complex formulation (an ex- 
ample for the CO2 molecule is presented in |27l1. How- 
ever, estimation of these vibrational parameters is rather 
difficult, the resulting values giving inaccurate results in 
some cases. In practice, a set of band-origin wavelengths 
and rotational constants for each vibrational transition is 
given, allowing the calculation of the line positions. 



The calculation of the lines intensities follows a similar 
approach than for diatomic spectra calculations. Emis- 
sion coefficients of rovibrational transitions of polyatomic 
spectra can be expressed through an expression similar to 
Eq. H 



47t" 



'SfijuFjiji 



1 1 v 



(15) 



Where the additional term Fjij/i designates the 
Herman- Wallis factor, which accounts for vibration- 
rotation interactions (see p. 110 in \'2ii\ ). This factor has 
been omitted from Eq. |3] as vibration-rotation interac- 
tions can be usually neglected in a rovibronic transition, 
owing to the usually large energy gap between the tran- 
sition electronic levels. However, for rovibrational tran- 
sitions, this interaction has to be accounted for, which 
explains it's inclusion in Eq. 1151 

The vibrational Einstein coefficient is related to the vi- 
brational transition moment squared $l v t v ti through the 
relation equivalent to Eq. 1121 for linear polyatomic rovi- 
brational transitions: 



64tt 4 
3hc* 1 



(2 -S ,y) 



(16) 



2.026 ■ 10~ 6 iA, 



(2-V) 
(2- V) 



Determination of this squared transition moment 
is rather complex. Instead, values of the integrated in- 
tensity of a vibrational band are tabulated at a refer- 
ence temperature To (usually 296 K) |29j . The value for 
the vibrational dipole moment (in atomic units squared 
(ean) 2 ) can then be deduced from the following expression 
1301137) : 



87T J 



■ S°, 



(ea ) 2 



y v'v" (2- * ,i")e*P 



hcE v , 



D 2 
(17) 

Expressions of the Honl-London factors for parallel and 
perpendicular transitions taken from |o II are presented in 
Tab. IO 

Finally, the Herman- Wallis coefficients can be devel- 
oped into the following polynomial expansions: 



P branch: (1 - Ax J" + A 2 J" 2 - A 3 J" 3 ) 2 
Q branch: (1 + A Q J"{J" + l)) 2 

R branch: (1 + A^J" + 1) + A 2 (J" + l) 2 + A 3 (J" + l) 3 
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A£ = A£ ^ 



where 



p (j"+e")(j"-e") (j"-i-£"Ae){j"-i"Ae) D=- — 

J" ..„ 2J" V 



(2,/" + l)l" 2 (,/" + !+£" Al)(./"-l"Al)(2,7" + l) 



J "(/" + 1 ) ,. ,, ^"^"t 1 ) „ V = i Az/ £ + J APf + 4AP2 



-n + + (J"+2+l"Al)(./" + l+l"Al) 2 

11 J" + l 2(J" + 1) 



Afj, 



Tab. 1: Honl-London Factors For Parallel and Perpendicular = 

Vibrational Transitions K 1 1<Ub5 + U - U47 — + 0.058-^ 



V 



C 2 = — p 

values for the different coefficients A being given for V (1.065 + 0.047 —y— + 0.058- 

each vibrational band. 



\' 2 



2.4 Line Shapes Calculations 

Line broadening mechanisms are grouped into two cate- 
gories: pressure broadening including natural, resonance, 
van der Waals, and collisional broadening mechanisms. 
Pressure broadening is represented by a Lorentzian pro- 
file for a given full width at half maximum (FWHM) A£>l, 
and Doppler broadening is represented by a Gaussian pro- 
file for a given FWHM Aug- A discussion on these dif- 
ferent broadening mechanisms can be found in 1321 . 

The convolution of the line shapes resulting from these 
two types of broadening mechanisms leads to a Voigt pro- 
file given by: 




However, the explicit convolution of these two line 
shapes can lead to high calculation times, and semi- 
empirical expressions approaching this exact line profile 
are seldom used. 

An empirical expression for this Voigt profile is given 
by Arnold [T3] 2 : 



v(y) = C\e 



+ 



1 + 4L> 2 



10 



W + D 2.25 

(19) 



2 the term 41n2 can be replaced by 2.772 for computational 
efficiency 



The accuracy to the exact expression being within 1 %. 

A simple modification to the Voigt FWHM expression 
has been proposed by Olivero lool: 



V = - ( 1.0692A^ L + 



0.86639A^£ 



4A*> 2 



giving a better accuracy of 0.02 % to the exact expres- 
sion. 

This presented lineshape calculation method presents 
the advantage of allowing a good accuracy to the exact 
function, without a great computational burden. How- 
ever, many other methods exist for the calculation of a 
Voigt lineshape. An overview and discussion for these 
different methods can be found in I34| . 

2.5 Strategies for Fast and Accurate 
Spectral Simulations 

Spectral calculations using the line-by-line approach 
can lead to large calculation times, as several thousands 
of lines are typically calculated and convoluted with the 
corresponding lineshapes. Also, such calculations can re- 
turn spectral grids which can be several million points 
wide. Therefore, some techniques may be used in order 
to reduce the number of calculated lines. A good way 
to achieve this consists in making some approximations 
regarding the fine-structure (spin-splitting effects) of the 
simulated spectra. This fine structure can in some cases 
be entirely or partially neglected resulting in lesser com- 
puted lines. 
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2.5.1 neglecting line spin-splitting effects 

Taking into account the effects of spin-splitting in cal- 
culations of the fine structure for multiplet transitions 
leads to the determination of up to 12 and 27 rotational 
branches for doublet and triplet transitions respectively 
instead of the 3 usual P, Q and R branches. This re- 
sults in calculation times increased by a factor of 4 and 
9 respectively. Therefore, one may want to neglect spin- 
splitting in order to reduce computation times, but with- 
out losing precision in the calculation result. One needs 
then to evaluate the importance of line spin-splitting over 
the line resolution (which depends on the lines FWHM). 
If the separation of multiplet lines is much smaller than 
their widths, such lines will be accurately modelled by 
one singlet line. 



2.5.2 Reducing the number of simulated rotational 
branches 

Even when spin-splitting effects need to be accounted 
for in spectral simulations, one may still reduce the num- 
ber of calculated rotational branches making use of cer- 
tain simplifications. 

The first case occurs for perpendicular transitions be- 
tween E and II states when neglecting the spin-splitting 
of the E state using the approach described in section 
12. 5. H leads to the superposition of several branches. 

As an example, for 2 E ^ 2 II doublet and for 3 E ^ 3 II 
triplet transitions 4 and 10 satellite branches respectively 
coalesce with the main branches when neglecting spin- 
splitting of the E state. 



This approach is usually valid for multiplet E states, 
which have small spin-splitting factors allowing transi- 
tions between such states to be accurately modelled as 
singlet states. 

A good example can be given for a 2 E — 2 E transition. 
For this case, the separation between two doublet lines 
(in A) is given by: 



AA ~ 10' 



(7«-7Gf)(J+g) 



(20) 



For the CN violet system , this is equivalent to a line 
splitting of 0.14 A for a rotational value as high as J=100. 
In classical spectrometry applications, this transition can 
be accurately simulated as a singlet transition, reducing 
calculation times without a loss of accuracy. For triplet E 
transitions, a slightly more complicated relationship can 
also be determined to evaluate whether the analyzed tran- 
sitions can be accurately modelled as singlet transitions 
(according to the transition 7 and A constants). 

The limits of this assumption occur when the spin- 
splitting of the studied states is no more negligible, and 
self-absorption of the spectra is high. In this case, a com- 
parison of a simulated spectrum using this simplification 
with a simulated spectrum accounting for spin-splitting 
shows non-negligible differences 1181 . However, even for 
multiplet states with larger spin-splitting values such as 
II and A, the same approach may still be valid for lower 
resolution spectra (with line widths of about 1-10 A or 
higher) . 

3 the spin-spin correction is 7 = 7.26 ■ 10 _3 cm~ 1 for the 
(X 2 E,v=0) state, and 7 = 17.16 ■ lO^cm- 1 for the (B 2 £,v=0) 
state 1351 



The second case occurs for calculations of high tem- 
perature spectra in applications such as combustion or 
plasma radiation in which some branches have small in- 
tensities and can be neglected. 

Energy exchanges between the different rotational 
states (R-R transfer) are very efficient which means that 
in most of the spectroscopic applications of interest, a 
Boltzmann distribution of these states may be assumed. 
Therefore, a line rotational intensity may be written as 



S{J) = 



(2 J + 1) exp 



F{J) 



qA' J' 
°A" J" 



Qrot 



2 J ■ 



(21) 



as each rotational state has a gj = 2 J + 1 degeneracy, 
the ground state will not be the most populated state. 
As an example, at room temperature, typically the most 
populated rotational level is around J — 7 — 10, whereas 
at a 2000 K temperature, the most populated level is 
around J = 20. 

This is an advantage for the simulation of high temper- 
ature gases radiation as the rotational transition prob- 
abilities (Honl-London factors) for some branches have 
smaller magnitudes than others, only reaching equivalent 
magnitudes for the first rotational levels (typically J < 5). 
If we take into account the fact that those levels are not 
the most populated according to what was said before, it 
can be verified that the peak rotational line strengths of 
such branches are typically more than one order of mag- 
nitude smaller than the peak rotational strengths of other 
branches. An example for a 3 E ^ 3 II transition at a 300 
K gas temperature is presented in Fig. [5] 

The ratio of the peak intensities between the second 
and first group of rotational branches is 3.3 for this case. 
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10 2 




Rotational Quantum Number 

Fig. 2: Rotational line strengths of the different branches of 
a 3 S <-> 3 fl transition for a Boltzmann distribution of 
the rotational levels at a characteristic temperature of 
300 K. Expressions for the Honl-London factors are 
taken from 1161 , The corresponding weak branches 
(group 1) are listed in appendix [B] 



For a 2000 K temperature, the same calculation results 
in a 15.8 ratio. Therefore omitting such weaker rota- 
tional branches for spectral calculations of high temper- 
ature gases has in most cases a negligible effect on the 
simulated spectra. 

Branches than can be neglected and coalesced branches 
when the spin-splitting of E states in perpendicular tran- 
sitions is neglected are presented in appendix [B] for dou- 
blet and triplet transitions between E and II states. 

3 Examples of Molecular Radiation 
Simulation in the Line-by-Line 
Approach 

Some applications of the different methods for the cal- 
culation of line positions are presented in this section. 
Specific spectra is numerically reproduced, and the appli- 
cability of each of these methods is discussed. The ex- 
perimental data presented here has been recorded at high 
resolution in low pressure plasma facilities, using emission 
spectroscopy techniques. 



3.1 Experimental Setup 

The molecular spectra presented in this study is issued 
from two similar low-pressure facilities where strongly 
emissive plasmas are obtained. The two facilities are lo- 
cated at the CORIA Laboratory in Rouen and the Labo- 
ratoire d'Aerothermique in Orleans respectively, both lab- 
oratories from the Centre National de la Recherche Sci- 
cntifiquc in France. Specific molecular systems were ac- 
quired using emission spectroscopy techniques, with typ- 
ical resolutions of 0.3-0.5 A. 



3.1.1 The SR5 Arc-jet Plasma Wind-tunnel 

The SR5 low pressure arc-jet plasma wind-tunnel avail- 
able at the Laboratoire d'Aerothermique has been used 
during the last years for the simulation of the entry con- 
ditions in Earth, Mars, and Titan planetary atmospheres 

A D.C. vortex-stabilized arc operating at low voltages 
(50-100 V) and low currents (50-150 A) delivers typical 
powers of 5-10 kW to the flow in the throat region of 
the nozzle. The low mass-flow rates (0.1-0.5 g/s) cross- 
ing the nozzle allow the obtention of an high-enthalpy 
(5-30 MJ/kg) steady plasma jet with a global yield of 50- 
70%. This level of specific enthalpy obtained with a small 
mass-flow rate is an advantage because of the low elec- 
trode erosion, which allows maintaining a steady plasma 
jet for several hours with a low level of contamination. 
The pumping system capacity of 26.000 m 3 h -1 ensures 
an ambient pressure of about 10 Pa to be maintained in a 
4.3 m long and 1.1 m diameter vacuum chamber in which 
the arc jet is expanded. 

Spectral measurements were carried using a SOPRA 
F1500 (Ebert-Fastie type) monochromator with a focal 
length of 1500 mm and a grating of 1800 grooves/mm 
sweeping a spectral region from 270 nm in the near- 
ultraviolet region to 950 nm in the near-infrared region. 
The grating is connected to an intensified optical multi- 
channel analyser (Princeton Instruments IRY 1024) . This 
device allows a 8.5-nm wavelength region to be expanded 
on 1024 pixels and is cooled by a Peltier element insuring 
an operating temperature of -35° C. The plasma is imaged 
onto the monochromator by a mirror telescope connected 
to the entrance slit by a quartz optical fiber. The en- 
trance slit opening can be adjusted, therefore modifying 
the experimental apparatus function which can reach val- 
ues down to a full width at half maximum (FWHM) of 
0.03 nm. 
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3.1.2 The CORIA Inductively Coupled Plasma Torch 

The CORIA inductively coupled plasma (ICP) torch 
is used for different research fields such as plasma de- 
pollution applications as well as the investigation of 
plasma radiation from the vacuum ultraviolet to the in- 
frared regions 1371 . 

The facility includes a high frequency (1.76 MHz) gen- 
erator, with typical voltages and intensities around 7 kV 
and 10 A respectively, producing an electromagnetic field 
onto a five fingers coil surrounding a quartz chamber in 
which a test gas is inserted. The typical mass flow rates 
of 1-5 g/s, for a delivered power of about 70 kW/s and 
a global yield of around 30 % lead to the formation of 
a high enthalpy (around 10 MJ/kg) stable and homoge- 
neous plasma. Further, a convergent nozzle can be added 
between the torch and aim long and 0.5 m diameter test- 
chamber in which measurements can be implemented. 

UV and visible spectra of air plasmas is acquired using 
an Acton Research spectrometer of 750 mm focal length 
implemented with a Princeton IMAX Intensified camera 
allowing spectral resolutions up to 0.045 nm. 

3.2 Presented Examples 

An investigation has been carried on the numerical re- 
production of molecular spectra from the CN Violet Sys- 
tem Aw = 0, the C2 Swan Bands Av = 0, and the 
First Negative System Av = 0. Each of these systems 
have been investigated previously 1261 138| and difficulties 
linked to the numerical description of such spectra have 
motivated their selection for a discussion in this work. 

3.2.1 The CN Violet System 

Conducting spectral measurements of line positions in 
low pressure and temperature gases and plasmas allows 
a precise and unambiguous line position determination, 
as broadening mechanisms (principally doppler broaden- 
ing, which is temperature dependent) will be kept to a 
strict minimum. However, compared to higher tempera- 
ture measurements, only the lower energy vibronic states 
will be excited, and less lines will be measured. There- 
fore, the interpolation of such line positions through the 
usual polynomial expressions will only be exactly valid 
for the lower vibrational and rotational levels. Extrapo- 
lation of such polynomial expansions to higher levels will 
of course be risky, specially taking into account the strong 
oscillations of the higher terms of such expansions. Con- 
sequently, a good balance between the accuracy of the 



line positions determination and the number of measured 
lines has to be achieved. 

The CN Violet system is a strong radiative emitter in 
most air and carbon species plasmas. Therefore, it has 
been extensively studied by different research teams for 
over 50 years, and various sets of spectroscopic constants 
have been proposed for the calculation of line positions. 
Such sets have been obtained through line position mea- 
surements to different levels of precision, and to different 
temperature ranges. Namely, fourier spectroscopy tech- 
niques have been applied to the measurement of this ra- 
diative system over the more recent years. 

Emission from the CN Violet system Av = has been 
measured in the SR5 facility for a low pressure and high 
enthalpy N2-CH4 plasma, in which a Boltzmann equilib- 
rium of the CN molecule rotational and vibrational levels 
is achieved. The apparatus function of the measurement 
system has been held at a FWHM of 1 A, so as to resolve 
solely the vibrational bands of the system and to allow 
an easier comparison of the different spectroscopic sets 
available regarding vibrational bandheads. 

The spectra has been numerically reproduced using the 
set of spectroscopic constants proposed by Herzberg |28| . 
the set of spectroscopic constants proposed by Huber & 
Herzberg )39l . the set of spectroscopic constants used in 
the 1985 version of the well known line-by-line spectral 
code NEQAIR 1401 . the sets of spectroscopic constants 
proposed by Cerny 03 for tne X 2 E + level and Ito 02] 
for the B 2 E + level, and finally the set of spectroscopic 
constants proposed by Prasad & Bernath 1351 . The tran- 
sition probabilities have been taken from the set proposed 
by Knowles 1431 . A comparison between the measured 
spectrum and the reproduced spectra, using the above 
constant sets and a line shape simulating the experimen- 
tal FWHM, is presented in Fig. El 

Large discrepancies can easily be detected between the 
different simulated spectra. All the spectroscopic con- 
stants reproduce to a good level of accuracy the 0-0, 1—1, 
and 2-2 bandheads of the experimental spectrum. How- 
ever, for higher vibrational bands, discrepancies tend to 
become larger, and some simulated spectra fails to re- 
produce accurately the band reversal observed in the ex- 
perimental spectrum which leads to the superposition of 
the 4-4 and 5-5 bandheads. Older spectroscopic constant 
sets such as the ones proposed by Herzberg and Huber & 
Herzberg are inadequate as they likely did not take into 
account higher vibrational levels when they were devel- 
oped, resulting in a limited validity range. 

The poor reproducibility of the measured spectrum by 
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Fig. 3: 



Comparison between a CN violet system Av=0 spec- 
trum measured in a CH4-N2 (0.002 g/s, 0.36 g/s, 
6.7 kW total power, 4 kW injected power, 3.8 Pa 
background pressure) plasma in the SR5 facility and 
simulated spectra using different spectroscopic con- 
stants. The molecule vibrational and rotational tem- 
peratures were iteratively set in the calculation until 
a best fit was achieved. The spectrometer apparatus 
function was determined using a mercury calibration 
lamp. The different parameters read: FWHM=lA, 
T„ = 10800 ± 100 K,T r = 4900 ± 100 K 



the simulated spectrum using the spectroscopic constants 
set of the NEQA1R85 code is likely due to the truncation 
to the value of oj e i/e for the vibrational levels polynomial 
expansion (see Eq. \'2:\t . Another example of such issues 



is discussed in [S] regarding the simulation of the N2 2nd 
Positive System. Care should therefore be exercised when 
using spectroscopic sets with large polynomial expansions 
in numerical codes which do not account for higher or- 
der terms. Resorting to the truncation of available sets 
may still yield correct approximations for the lower vibra- 
tional levels energies, but for higher levels energies, which 
are more dependant on higher order terms of such poly- 
nomial expansions, this is no longer the case. Instead, 
it is preferable to fit the initial polynomial expansion to 
a lower order polynomial expansion, using the validity 
range of the initial fit for the level energies calculations. 

More recent spectroscopic datasets, determined 
through spectral measurements using Fourier spec- 
troscopy techniques, reproduce more accurately the ex- 
perimental spectrum. The spectroscopic constants pro- 
posed by Cernyi and Ito already allow reproducing cor- 
rectly the 4-4 and 5-5 bandhead superposition. How- 
ever, the 6-6 bandhead still does not match accurately 
the experimental spectrum. Finally, the most recent spec- 
troscopic constants, proposed by Prasad allow an exact 
reproduction of the experimental spectrum (to the mea- 
sured resolution). 



3.2.2 The C2 Swan Bands 

Various sets of equilibrium spectroscopic constants are 
proposed by different authors for the simulation of differ- 
ent diatomic spectra. For most applications, the accuracy 
of such datasets is sufficient for carrying spectral simula- 
tions, provided that an analysis of the available data is 
carried and the more accurate data is selected. However, 
in some minor cases equilibrium spectroscopic constants 
no longer suffice for a correct description of measured 
spectra. The frequently observed Swan Bands of the C2 
molecule provide a good example for when equilibrium 
spectroscopic constants fail to accurately reproduce ex- 
perimental spectra. 

A good starting point for the analysis of the spectro- 
scopic constants for the upper and lower electronic states 
of this system is given by a study from Phillips |44|. which 
analyzed the first 10 vibrational levels of such states. Al- 
though the vibrational levels of the a 3 n u state follow the 
usual monotonous evolution pattern typical of Eqs. I2al 
I2bl and l2cl the vibrational levels of the d 3 n 9 state show a 
consistent change of trend, starting from the vibrational 
level v=5. In this specific case, it is no longer possible 
to fit the level spectroscopic constants G v , B v , D v , ... to 
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Vibrational Quantum Number 



of Eq. US is no longer adapted to the description of the 
level-dependent rotational constants B v , yielding a large 
fit dispersion. 

The origin of the perturbation observed for the spec- 
troscopic constants of the d 3 II 9 electronic state, starting 
from level v=5 has been investigated through an analysis 
of the reconstructed potential curves of the transition and 
the neighboring triplet states, using the RKR method. 
Spectroscopic constants for the b 3 E 9 and e 3 n g states 
were taken from the compilation of Huber & Herzberg 
|39| . and the spectroscopic constants for the c 3 E+ state 
were taken from Ballik and Ramsay 1491 who determined 
them through the study of perturbations in the A 1 U U 
state. The calculated potential curves are presented in 
Fig. 



Comparison between level and interpolated rotational 
constants B v for the a?Tl u and the d 3 II 9 electronic 
states of C2. A vibrational perturbation of the d 3 II 9 
state leads to a trend change for the rotational con- 
stant above the v=5 level, preventing an accurate fit 
to the usual equilibrium spectroscopic constants from 
being obtained 




Fig. 5: 



Internuclear Distance (A) 



Potential curves of the d 3 n 9 , a 3 Il„ and neighbor- 
ing triplet electronic states obtained using the RKR 
method. The dashed lines represent the shape of the 
d 3 n 3 and e 3 II 9 states potential curves if the vibra- 
tional perturbation resulting from the avoided inter- 
section of the potential curves is not accounted for 



the usual polynomial expansions. A fit of the rotational 
constants B v provided by Phillips to the fifth order poly- 
nomial expansion of Eq. I2bl has been carried for both the 
electronic levels, and the interpolated values for B v are 
compared to the initial values in Fig. 0] 

It can be easily verified that, although the values of B v 
for the a 3 II u electronic state can be exactly reproduced by 
a polynomial expansion, this is no longer the case for the 
d 3 II s electronic state where the mathematical expression 



An analysis of the curve plots show that the perturba- 
tion of the d 3 II 3 electronic state could be explained by 
the crossing of the b 3 £ 9 potential curve which occurs at 
v=5. However for such cases, the perturbation would re- 
main local, and higher vibrational levels of the d 3 IT 9 state 
would remain unperturbed, which is not the case. There- 
fore, the more probable cause for this perturbation could 
result from an avoided intersection between the d 3 II 9 and 
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e 3 II g potential curves (see p. 295 in |28p. which alters the 
curve shape above the perturbation. 

Dashed curve lines have been added to the plot to show 
the avoided intersection between the potential curves. 
The reconstruction of the unperturbed potential curve 
for the d 3 n a electronic state has proceeded including 
solely the unperturbed levels t>=0-4 in the RKR calcu- 
lation. The extrapolation of the potential curve by a 
Hulburth and Hirschfilder expression has taken into ac- 
count a higher dissociation level [C( 3 P)+C( 1 D)]. The re- 
construction of the e 3 n g electronic state potential curve 
has proceeded identically, taking into account a lower dis- 
sociation level [C( 3 P)+C( 3 P)]. 



An experimental spectrum from the C2 Swan Bands 
Av=0, obtained in the CORIA ICP torch for a low pres- 
sure C02-N2-Ar plasma flow, has been simulated using 
the level constants proposed by Phillips, and their inter- 
polation to equilibrium constants. The transition proba- 
bilities used in the calculation are issued from Cooper 1 151. 
The comparison of the simulated spectra, using level and 
equilibrium spectroscopic constants to the experimental 
spectrum is presented in Fig. |5] 

As it could be expected after the analysis of the fits 
presented in Fig. 01 it is verified that equilibrium con- 
stants, although fitted to very accurate level constants, 
fail to reproduce exactly the experimental spectrum. This 
is particularly apparent for the 1-1 bandhead at 5120 A, 
which is not very well reproduced unlike the spectrum is- 
sued from level spectroscopic constants. Therefore, when- 
ever vibrational perturbations are present for any of the 
states of the radiative transition, it is likely that the ex- 
plicit insertion of level spectroscopic constants in simu- 
lation codes will be a required feature for a sufficiently 
accurate simulation. 



- Simulation With Equilibrium Constants 



3.2.3 



The Nt First Negative System 



A final example can be presented for a spectrum af- 
fected by rotational perturbations of it's internal states 
(intersection of potential curves for J / 0). For this case, 
there is a shift of the nearby levels populations following 
a mathematical expression of the type 1/x (see p. 283 in 
I28| l. The more accurate but more complex method for 
accounting for the presence of such perturbations consists 
in solving the perturbative hamiltonian for the system 
However, in our case, it was chosen to resort to the usual 
Klein-Dunham coefficients for the calculation of the level 



AAAAAA 





5140 5150 
Wavelength (A) 



Fig. 6: Comparison between a C2 Swan bands Av=0 spec- 
trum measured in an CO2-N2— Ar (2 g/s, 0.2 g/s, 
0.2 g/s, 70 kW total power, 20 kW injected power, 
1000 Pa background pressure, no nozzle) plasma in the 
CORIA ICP torch using equilibrium and level spectro- 
scopic constants. The molecule vibrational and rota- 
tional temperatures were iteratively set in the calcu- 
lation until a best fit was achieved. The spectrome- 
ter apparatus function was determined using a mer- 
cury calibration lamp. The different parameters read: 
FWHM=0.55 A, T„ =T r = 4200 ± 100 K 



energies, with the addition of a term in the form of 



Ej = Ej 



AE„ 



2 (J — Jpert — 1/2) 



(22) 



The upper electronic state B 2 £+ of the First Neg- 
ative System is known to be perturbed by the A 2 n u elec- 
tronic state. The influence of the shifts of levels energies 
on the observed spectrum are quickly appreciable, which 
makes this radiative system a good candidate for the il- 
lustration of the importance of perturbations modelling 
in spectral simulations. Such perturbations have been 
extensively studied in |47l I46| . and the values for J pe rt 
and AE ra „ (position and intensity of the perturbation) 
for the perturbations resulting from the above described 
interaction are reported in the references therein. 

Emission from the First Negative System system 
Av = has been measured in the SR5 facility for a low 
pressure air plasma. The apparatus function of the mea- 
surement system has been set at a FWHM of 0.3 A, which 
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allows fully resolving the rotational lines of the spec- 
trum. Spectral simulations have then been carried for 
this molecular system, firstly excluding, and lately tak- 
ing into account the influence of the perturbations in the 
spectral model. Level spectroscopic constants provided 
by Michaud I47| and transition probabilities provided by 
Laux |48| have been taken into account for the calcula- 
tions. The obtained results and a comparison with the 
experimental spectrum are presented in Fig. |7| 



- Simulation With Deperturbed Constants 




- Simulation With Perturbed Constants 
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Fig. 7: 



Comparison between a 1st negative system Av=0 
spectrum measured in an air (0.2 g/s, 4.6 kW total 
power, 2.4 kW injected power, 4.4 Pa background 
pressure) plasma in the SR5 facility and simulated 
spectra with and without the simulation of perturba- 
tions. The molecule vibrational and rotational tem- 
peratures were iteratively set in the calculation un- 
til a best fit was achieved. The spectrometer ap- 
paratus function was determined using a mercury 
calibration lamp. The different parameters read: 
FWHM=0.33 A, T„ = 4000 ± 100 K, T r = 2600 ± 100 K 



The analysis of the simulated spectra allows detecting 
large discrepancies due to the presence of perturbations. 
Namely, the structure of the 0-0 vibrational band is com- 
pletely altered by the perturbations of the upper levels 
energies. The effects of the perturbations in the simu- 
lated spectrum are very apparent near 3900 A. It can be 
verified that the simplified model used for including the 
effects of perturbations in the level energies allows ob- 
taining a very close match to the experimental spectrum. 
This is no longer the case for the unperturbed spectrum, 



which fails to reproduce correctly the experimental spec- 
trum, although accurate level spectroscopic constants are 
used in the calculation. 



4 Conclusion 

The methods and issues regarding line-by-line spectral 
calculations have been discussed in this paper. Accurate 
methods for calculating the main parameters of bound 
radiation such as line positions, intensities, and shapes, 
have been presented for diatomic rovibronic transitions, 
and linear polyatomic rovibratinal transitions. 

Moreover, as line-by-line molecular simulations lead to 
the calculation of spectrums composed from a very large 
number of lines, with arising issues regarding computa- 
tion times and memory limitations, strategies leading to 
a reduction in the number of calculated lines, with minor 
accuracy losses, have been presented. 

A more complete discussion on the different methods 
available for the calculation of line positions has also 
been carried. Additionally to the necessity of selecting 
the most adapted method for the calculation of line po- 
sitions, the importance of selecting the most accurate 
and adapted spectral dataset has been highlighted. It 
has been showed that coping with the specific issues of 
the simulation of unperturbed but also perturbed spectra 
may allow reaching simulated spectra that matches very 
accurately experimentally determined spectra. Moreover, 
achieving such accuracies in spectral simulations also re- 
sults in a very suitable method for determining vibra- 
tional and rotational temperatures T v and T r to a good 
level of accuracy. The methods which were applied to the 
simulation of the presented experimental spectra have al- 
lowed determining the vibrational and rotational temper- 
atures to an accuracy of 100 K. Even when a Boltzmann 
equilibrium is not achieved for such levels, such methods 
can be successfully applied to the determination of the 
levels population distributions |53| . 
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Appendix 

A Multiplet Level Energies 

All the expressions below are written in wavenumber 
units (cm -1 ). 

Doublet Levels 

The general expression for the doublet energy levels is 
given by |28| : 



IT states energy levels are given by 



j n 2 (J>2) = B v 



D v J 



J (J + 1) - y/ yi +4J(J + l) 

2 y 2 -2J(.J+l) 

3 hi+4J(J+1) 

4 



1 



(26a) 



IIi(J>l) =B V 



J(J + 1) + 



4 y 2 - 2 J( J + 1) 



3yi+4J(J + l) 



F 3 /2{J>1) =B V 

-D V J 
F 1/2 [J>0) = B v 



(J+|) 2 -A 2 

1\2 



4 ( J +§) +Y(Y-A)A< 



2\ 2 



D V [J + 



(26b) 



TI (J>0) = B V 



(23a) 



(J+|) 2 -A 2 

1\2 



+1 U{J+±) +Y(Y-4)A< 



2\ 2 



J(J+1) + VV+4J(J+1) 

2 j/ 2 -2J(J+l) 

3 hi+4J(./+1) 

4 



A, ( J + 5 



(26c) 



A,(J + 1) 4 



(23b) 



For E states, a simpler expression of the levels energies 



is used EH]: 

2 £ 3 / 2 (J>l) = B U (J(J + 1)) - D V {J(J+ l)) 2 + 7 
2 E 1/2 (J>0) = B„(J(J + 1)) - D V (J(J + l)) 2 - 7 



with 



2/i = Y(Y-4) + - 



V2 =Y(Y-1)-- 



r 

(24a) 

j + 1 \B Coalesced, Weak and Strong 

2 J Branches for Doublet and Triplet 

(24b) Transitions 



y _ ^ 

B v 



Triplet Levels Expressions for parallel 2 E ^-> 2 E transitions for Hiind 

The more accurate expressions for 3 E states energy lev- case b are taken from G3> expressions for parallel 3 E • 

no _ E transitions for Hiind case b are taken 

| are preferred to the ones 

found in 



els given by [501 form the formulae of EH with a typo- E transitions for Hiind case b are taken from HJ, expres- 
graphical correction from are preferred to the ones sions for perpendicular 2 E ^ 2 IT and 3 E ^> 3 n transi- 
tions for the intermediary case between Hiind cases a and 



j E 2 (J>2) = B V {J(J + 1)) - D V {J{J + 1)Y 
- ( \ v - B v + i 7u ) 

, 2 



b were taken from 1181 and I16| respectively. Expressions 
for parallel 2 n <-> 2 n and 3 n <-> 3 n transitions for the in- 
(25a) termediary case between Hiind cases a and b were taken 
from HFJ. 



[ (A„ - B v + i 7 „) 2 + 4 J( J + 1) (B„ - i>) 2 



s Ei(J> 1) = B V {J(J + 1)) - D V {J{J + l)) 2 
3 E o (J>0) = B V (J(J + 1)) - D„(J(J + l)) 2 
- ( \ v - B v + i 7u ) 

, 2 



(25b) Weak Rotational Branches 



(25c) 



For multiplet parallel transitions (AA = 0), all the 
satellite rotational branches have weak strengths. For 
-|n|ultiplet perpendicular transitions (AA = ±1) the weak 



(X v -B v + |>) + 4J(J + 1) [B v - i>) Jrotational branches are listed in Tab. H 
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Q P21 


R P 3 1 


°Pl2 


Q Rl2 


R Q21 
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R 3 1 


Q Rl2 


S i?21 
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Q P 3 2 


N Pis 




P Q23 


R Q32 


°Ql3 




Q i?23 


S i?32 


R13 



Tab. 2: Weak rotational branches for doublet and triplet 
transitions 



List of Coalesced Branches for E — II 
Transitions Without Spin-Splitting Effects 
for the S State 

The coalesced rotational branches for multiplet E <-» II 
transitions with a neglected spin-splitting of the E state 
are listed in Tab. El 



l E - 2 n 



m - 2 s 



n 



P 0l2 = p 2 

Q P 2 i = Qi 
Q Ri2 = Q 2 
R Q2i = Ri 



P Q 



12 



Q Rl2 
Q P21 
R Q21 



Pi 

Qi 
Q2 
R2 



R31 



R Q 



S R21 
°Pl2 
?32 : 
3 23 
P Ql2 
N Pl 3 - 

P Q 23 



S R 3 



Pi 
Pi 

P2 

Qi 
Q2 
Q2 
Q3 
R2 
R 3 
R3 



R Q21 ■■ 
Rt 



32 
31 
12 
23 
21 
•32 
N Pl3 
P Ql2 



°Pl 
°P, 

S R 



Pz 

= P 3 
= P 3 

= Qi 
= Q 2 
= Q 2 

= 0.3 

= Ri 
= Ri 

= R2 



Tab. 3: Coalesced rotational branches for multiplet E — II 
transitions when neglecting spin-splitting for the 
E state 
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